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ABSTRACT 

We study the gravitational lensing effects of spiral galaxies by taking a model of the Milky 
Way and computing its lensing properties. The model is composed of a spherical Hernquist 
bulge, a Miyamoto-Nagai disc and an isothermal halo. As a strong lens, a spiral galaxy like 
the Milky Way can give rise to four different imaging geometries. They are (i) three images 
on one side of the galaxy centre ('disc triplets'), (ii) three images with one close to the centre 
('core triplets'), (iii) five images and (iv) seven images. Neglecting magnification bias, we 
show that the core triplets, disc triplets and fivefold imaging are roughly equally likely. Even 
though our models contain edge-on discs, their image multiplicities are not dominated by disc 
triplets. The halo is included for completeness, but it has a small effect on the caustic structure, 
the time delays and brightnesses of the images. 

The Milky Way model has a maximum disc (i.e., the halo is not dynamically important 
in the inner parts). Strong lensing by nearly edge-on disc galaxies breaks the degeneracy 
between the relative contribution of the disc and halo to the overall rotation curve. If a spiral 
galaxy has a sub-maximum disc, then the astroid caustic shrinks dramatically in size, whilst 
the radial caustic shrinks more modestly. This causes changes in the relative likelihood of the 
image geometries, specifically (i) core triplets are now ^9/2 times more likely than disc 
triplets, (ii) the cross section for threefold imaging is reduced by a factor of ~ 2/3, whilst 
(iii) the cross section for fivefold imaging is reduced by ~ 1/2. Although multiple imaging is 
less likely (the cross sections are smaller), the average total magnification is greater. The time 
delays are smaller, as the total projected lensing mass is reduced. 

Key words: gravitational lensing - Galaxy: structure - Galaxy: disc - Galaxy: bulge - 
Galaxy: halo 



1 INTRODUCTION 
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stratified on similar 



In gravitational lensing, galaxies 
very idealized manner. For example 
the potential or the density 

trie ellipses have been widely studied iKassiola & KovneJI 
Kormann. Schneider & Bartelmannlll994IWittlll996l : IWitt < 
19971 : iHunter & Evansll200lhlEvans & HunterlbOtm Such simple, 
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analytic models often yield valuable insights and are reasonable 
enough as representations of elliptical galaxies. They are less well- 
suited for spiral galaxies, in which the effects of the disc, bulge and 
halo all need to be taken into account. 

Gravit ational lensing by spiral g alaxies has received atten- 
tion from iKeeton & Kochanekl ( 1199 8l) . who studied the proper- 
ties of Mestel and Kuz min discs embedded in isothermal haloes. 
IWang & Turned Jl997l) found an analytically tractable model of 
a homogeneous disc embedded in an isothermal sphere, which is 
nonetheless realistic enough to illustrate a number of the important 
effects. The presence of a disc makes lensing effects sensitive to 
inclination. In particular, a nearly edge-on disc causes three images 
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to form on the same side of the lens, a config uration which ha s 
possibly been observe d in APM08279+5255 dlbata et alj[l999h . 
Moll er & Blainl { 19981) used ray-tracing to con firm some of the the- 
oretical results of Keeton & Kochanek (1998). 



Of the ~ 90 stro ng lenses, only a handful are known to be spi- 
ral galaxies (see e.g.. Ilackson et al]|200d : Iwinn et alj|2003l) . The 
most ext ensively studied spiral galaxy strong lens is B 1600+434 
(see e.g., Koopmans 'et alJl998l : |jackson et alfe oOO). This is a two- 
image system, composed of a nearly e dge-on lens at z = 0.414 
lensing a background quasar at z = 1.6. IWinn et alj f2003h studied 
another two-image system, PMN J2004-1349, for which the lens is 
an isolated spiral galaxy with an inclination angle of ~ 70°. There 
are at least two known four-image lenses, Q 2237+0305 (see e.g., 
ISchmidt et ai1ll998l: iTrott & Webstej|2002l) and B0712+472 (see 
e.g., Jackson et alJl998l ; lKawano et al.ll2004l) . lBartelmann & Loebl 
( 1998) suggest that the underabundance of spiral lenses in optical 
searches may be caused by extinction by dus t and that the miss- 
ing lenses can be recovered in radio searches. Bartel manr] J200fj|) 
makes the important point that lensing by spiral galaxies may dom- 
inate over ellipticals once angular resolutions of the order of 0.1" 
can be achieved. He showed that for the projected Next Generation 
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Space Telescope, there may be of the order of 10 quasars per square 
degre e brighter th a n V ~ 26 l ensed by spirals. 

iMaller et al.1 feOOOh and IWinn et al.l J2003h emphasise the 
unique contribution that studies of gravitational lensing of edge-on 
spirals can provide. In the Milky Way galaxy, it has been a mat- 
ter of debate as to whether the disc is "max imum" or not (see 
e.g.. lSellwood & Sanderll988l ; lsltckettlll997l) . A spiral galaxy has 
a maximum disc if the inner parts of its rotation are completely 
dominated by the disc and bulge, rather than the dark halo. Al- 
though recent evidence strongly supports the idea o f a maximum 
disc for the Milky Way (Englm aier & Gerharcjl 19991 ; Ifflfner et al.l 
2000; Binnev & Evang l200ll) . the relative contributions of the 
disc, bulge and the halo to the flat rotat ion curves of external 
spiral galaxies remains un certain (see e.g. ICourte au & Rixj |l999l : 
Palun as & Williams! [2000). The issue of whether spiral galaxies 
typically possess a maximum disc could be resolved by the study of 
a sample of edge-on spirals acting as lenses, as gravitational lens- 
ing effects are sensitive to the relative masses in the rounder halo or 
flatter di sc componen t s. The first steps in this direction have been 
taken bv lMalleretaII ( l2000h who concluded that the spiral lens of 
B 1600+434 is not maximum. 

In this paper, we pick a model of the Milky Way galaxy 
with bulge, dis c and halo that is wi d ely-used in Galactic astron- 
omy (see e.g., Johnston et alj 1 19951 ; iDinescu et al.l 1 19991 ; iHelmil 
|2004| ; iFeUhauer et alj|2006t) and ask what its lensing properties 
would be, if it were acting as a gravitational lens at a redshift of 
z ~ 0.4. The model is composed of a Hernquist (1990) bulge, 
Miyamoto-Nagai (1975) disc and an isothermal halo, combined so 
as to reproduce a flatfish rotation curve of amplitude ~ 220 kms" 1 
at the Sun. This model has a maximum disc. We first discuss the 
lensing properties of the Hernquist bulge and Miyamoto-Nagai disc 
in §2 and §3, before combining them with an isothermal halo in §4 
to give the full Milky Way model. We investigate how the lensing 
properties change as the masses of the components are varied - in 
particular, as the disc changes from maximum to sub-maximum. 



2 THE HERNQUIST BULGE 

The Hernquist (1990) model was developed for elliptical galaxies 
and the bulges of spiral galaxies. The three-dimensional Hernquist 
mass distribution is 



p(f) 



2-7T f(f + ro) 3 



(1) 



where ro is a "core radius" and f is the spherical polar radius. Inte- 
grating l[T|l along the line of sight yields 



2tt r 2 (l-u 2 )< 



[(2 + u 2 ) X (u)~3,] , 



(2) 



where R is the radial coordinate in the lens plane and u = R/ro, 
while 
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For convenience, we follow Schneid er et al. ||l992|) in defining 
the dimensionless quantities 
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Figure 1. Critical curves (left) and caustics (right) of Hernquist lenses with 
different strengths of shear. The regions of multiple imaging in the source 
plane are labelled. In the top panels (7 P = 0), spherical symmetry ensures 
that the tangential caustic is a degenerate point at the origin. In the middle 
panels (7 P = 0.1), the tangential caustic (full curve) is wholly contained 
within the radial caustic (dotted curve). In the lower panels (7 P = 0.2) , the 
cusps of the tangential caustic pierce the radial caustic (so called "naked 
cusps"). 



with D s , D\, D\ s being the distances to the source, lens, and be- 
tween lens and source, and £o an arbitrary scale length. Now, the 
bending angle for a circularly symmetric lens is 



, , 2 f r , , , m(r) 
a(r) = — I r K(r ) ar = 



(5) 

Integrating (|2} and choosing the scaling fo = r o gives us the bend- 
ing angle of the Hernquist model as 



a{r) = m b 



r(xM-l), 



m b 



(6) 



Introducing k(t) = m(r) /r 2 , then the lens equation including any 
quadrupole perturbation is 



= [1 - R(r)] 



Ti 






r 2 



(7) 



where (x,y) are dimensionless (length scale £o) Cartesian coor- 
dinates in the lens plane and (£, rj) are dimensionless coordinates 
in the source plane (length scale r/o = fo-D s / ' D\). The compo- 
nents ri,2 are given by k p ± 7 P , where k p is the density of a uni- 
form sheet of matter superimposed on the lens and 7 P is the shear. 
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Figure 2. Level curves of the logarithm of the convergence (log 10 k) for the 
Miyamoto-Nagai disc with raj = 1.0. Contours are at intervals of 0.5. The 
projected density can be dimpled (upper panels) or highly flattened (lower 
panels). 
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Figure 3. The central ellipticity of the projected mass distribution of the 
Miyamoto-Nagai disc is constant along lines of constant gradient in (a, b) 
parameter space. Along such lines, models closer to the origin are 'diskier' 
(possess a stronger ridge along the x-axis). 
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Direct differentiation yields the components of the Jacobian ma- 
trix Aij = t^-, from which we obtain the image magnification 

^ = 1 / det A , 



An = I-kM-Ti- — k'M 

r 

XV- i 



An = An=-— R'(r) 
r 



(8) 



A 22 = l-R(r)-r a -^R'(r). 

r 

The critical curves are given by /i" 1 = det A = 0, i.e. 

2 2 

(i-s-ri)(i-«-r 2 )-(i-/t-ri)^-/t'-(i-K-r 2 )— «' = o.(9) 

r r 

Neither the critical curves nor the image positions can be found an- 
alytically, but both can be found numerically. The critical curves 
and caustics are plotted in Figure [T] for different 7 P for mj, = 
1, ftp =0; the numbers of images of sources in different regions of 
the source plane are shown. Varying mh and k p changes the sizes 
of the critical curves and caustics but not their asymmetry. The 
number of images is odd, not even. The three-dimensional mass- 
density is singular at f = and the surface density is logarithmi- 
cally divergent at r = 0, but the bending angle a and the deflec- 
tion potential are continuous and finite, so the od d-number theorem 
applie s. This is consistent with the study of lEvans & Wilkinson! 
dl998h on cusped density distributions and image numbers. They 
found that that cusps less severe than isothermal oc r _1 ) give 
rise to an odd number of images, but those more severe than isother- 
mal give rise to an even number of images. 
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Figure 4. Critical curves (left) and caustics (right) for the Miyamoto-Nagai 
disc with raj = 1.0 and various values of a and b. The central ellipticiy 
of the projected density contours is 0.16 (upper panels), 0.37 (middle) and 
0.67 (lower). The regions of multiple imaging in the source plane are la- 
belled. Note that the disc model contains the same mass, but that the central 
surface density is diminishing on moving from the top panels to the bottom 
panels. 
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Figure 5. Rotation curve for the model of the Milky Way built with disc, 
bulge and halo (solid line). The dashed line shows the effect of removing 
the halo component, whilst the dotted line shows the sub-maximum model. 



3 THE MIYAMOTO-NAGAI DISC MODEL 
3.1 Preliminaries 

The Miyamoto-Nagai (1975) disc has the following density distri- 
bution: 



4n[R 2 +(A+VB 2 +Z 2 ) f/ 2 [B 2 + Z 2 ] 3 / 2 



Here, R and z are cylindrical coordinates. The A — *• and B — *• 
limits are 



p(R,z) = 
p(R,z) = 



35 Ma 



4tt (R2 + z 2 + #2)5/2 
AM d 1 



^ (R 2 + A 2 ) 



3/2 



S(z). 



(11) 
(12) 



which are the Plummer (1911) model and the Kuzmin (1956) 
disc. The lensing prop erties of the form er are discussed in 
IWerner & Evansl d2006l) and the latter in iKeeton & Kochanekl 
dl998h . 

The Newtonian potential of the Miyamoto-Nagai disc is 



-GM 



R 2 + [A+ V / z 2 +B* 



-1/2 



(13) 



which is obtained from that of a point mass by the substitution 

2 



f 2 = R 2 + z 2 



R 2 + [A+yJz 2 + B' 



(14) 



3.2 The Projected Matter Density 

The integral of the Newtonian potential along the line of sight 
(which is proportional to the deflection potential) and the integral 
of the 3-D density along the line of sight (which is the surface mass 
density relevant to lensing) are not generally analytic. However, in 
the most important case when the disc is viewed edge-on, these 
integrals are analytic. From 

/+oo 
p(x',y',z')dy' (15) 
-OO 

we obtain the dimensionless surface mass density 
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Figure 6. Dependence of <r cr on lens and source redshifts in an Einstein-de 
Sitter universe. Contours are in units 1O 9 M0 kpc -2 . 



_ b 2 m d [2y/b 2 +y 2 (2a 2 + b 2 + y 2 ) + a(a 2 + 5b 2 + x 2 + 5y 2 )] 
2(b 2 + y 2 )' i / 2 [a 2 + b 2 +x 2 +y 2 + 2a^b 2 + y 2 ] 2 
where we have set 

M d 



(16) 



x z 
so to 



m d = 



(17) 



i^2 ' 

A h B 

so so 

so that (x,y) are dimensionless Cartesian coordinates in the lens 
plane. 

Some level curves of log k are shown in Figure|2] The model 
is useful, as the surface density distribution can be highly flattened. 
We see that a controls the overall extent of x-y asymmetry of the 
surface mass density (a = is circularly symmetric), while b con- 
trols the sharpness of the ridge along the x-axis. Asymptotically, k 
falls off along the x- and y-axes as 

+ 0(x" 4 ) , 



, m d a _ 2 

K(X, 0) = — X 



2b 



K(0,y) = b 2 m dy - 4 + O(y- 5 ) 



(18a) 
(18b) 



Near the centre, the level curves of k are ellipses with semiaxes 
| O20 1 and | ao2 1 , where 
a + 4b 



a,20 = —m d 



(Hv> 



-ma 



26(a + 6) 4 ' 
3a 2 + 9ab + i 



and 



4b 3 (a + b) 3 

are coefficients in the Taylor expansion about the origin 

k(x, y) = k(Q, 0) + a 2 QX 2 + a 02 y 2 + 0(x 4 , y 4 ,x 2 y 2 ) 



(19a) 
(19b) 

(20) 



The isodensity contours (or isophotes given a constant mass-to- 
light ratio) of the model near the origin are the ellipses given by 
l !19t and < |20b - The central ellipticity is 



2(c + 4) 



(c+ l)(3c 2 +9c + ! 



(21) 



where c = a/b. This vanishes when c = 0, as it should. Contours 
of constant central ellipticity are straight lines in the (a, b) plane, as 
shown in Figure [3] Moving along a line of given central ellipticity 
eo away from the origin takes us from models that are disc-like in 
the outer parts to models that are intrinsically rounder. 
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3.3 Disc lensing 

The deflection potential, related to the surface mass density by 

2 



K ( x ,y) = |v 2 i/> , is 



x 2 + ( a + vV + b 2 



V> = - m d log 



Note that the b — » limit is 

^ = \ mA lo s [ x ' 2 + ( a + \y\) 2 ] > 



(22) 



(23) 



not to be confused with the deflection potential of a point mass at 
y = —a, which is identical except \y\ — > We remark that the 
very simple projected potential of the Miyamoto-Nagai disc seems 
not to have been noticed before. It is an attractive model for studies 
of lensing of arbitrarily flattened mass distributions. 

The bending angle a — Vi/> is algebraic, with components 

ma x 



x 2 + 



(a + y/y 2 + ft 2 ) 

(a + y/y 2 + b 2 



m d y 



vV + b 2 ( x 2 + (a + y/y 2 + b 2 



(24) 



The critical curves and caustics can again be found numerically, 
solving 



detA = (1 - xjj xx ){l - xjj y y) - {-i> X y) 2 







(25) 



by Newton-Raphson along radial lines. A sample showing qualita- 
tively different configurations is shown in Figure [4] with the num- 
ber of images for sources in different regions marked. In the upper 
panel (eo = 0.16), the tangential caustic is contained within the 
radial caustic. In the middle panel (eo = 0.37), the cusps of the 
tangential caustic are "naked". The density contours become flatter 
still in the lower panel (eo = 0.67), and the tangential caustic has 
vanished. This is easily understood, as although the disc mass is the 
same, the central density diminishes on moving from the top to the 
bottom panels. Hence, this i s the sequence from the tw o-lips to the 
single-lips caustic (see e.g jKassiola & KovneJ[l993l) . 



Table 1. Mean magnifications and time delays (in days) for the core triplet, 
disc triplet and five image geometries for the Milky Way Model, the sub- 
maximum disc, and the Milky Way model with halo removed 





Milky 


Sub-Maximum 


Halo 




Way 


Disc 


Removed 


5-image systems 








Total magnification 


4.8 


7.2 


3.6 




2.5 


3.6 


1.9 




1.8 


2.7 


1.4 


^3 


0.37 


0.72 


0.30 




0.14 


0.26 


0.12 


Ml 


0.005 


0.007 


0.006 


*5-tl 


39.0 


23.9 


35.4 


*4 — t\ 


31.6 


18.8 


29.1 


t 3 -tl 


10.8 


3.9 


11.4 


t 2 ~tl 


8.5 


3.0 


8.7 


Core triplets 








Total magnification 


2.1 


2.6 


1.7 


M3 


1.6 


2.0 


1.2 




0.51 


0.59 


0.42 


Ml 


0.02 


0.02 


0.02 


tz-H 


56.8 


34.8 


51.4 


ti-t\ 


51.3 


31.5 


46.6 


Disc triplets 








Total magnification 


4.2 


6.2 


3.4 


M3 


1.9 


2.7 


1.5 


M2 


1.6 


2.3 


1.3 


Ml 


0.69 


1.17 


0.58 


tz-ti 


1.27 


0.30 


1.54 


t2-h 


1.04 


0.24 


1.30 



4 THE MILKY WAY AS A STRONG LENS 
4.1 A Model for The Milky Way 

A widely used model of the Milky Way is a combination of a 
Hernquist bul ge, a Miyamoto-Nagai disc and a cored isothermal 
halo (see e.g. IPaczvnskil 19 90; Johnst on et alJl995l ; lDinescu et al.l 
1999). The Newtonian potential (up to an additive constant) is 



$ = $ b + $d + Oh 



(26) 
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Figure 8. Fermat surfaces and image positions for the sources marked 1 to 4 
in Figure|7] Top left: Source 1 produces a 'disc' three-image configuration, 
with two brighter images straddling a demagnified one in the plane of the 
disc. Top right: Source 2 produces a three-image geometry. The central im- 
age is highly demagnified. Bottom left: Five-image geometry from source 3. 
Again, the images in the plane of the disc are demagnified, with the central 
one highly demagnified. The images at (4,4.6) and (5,-2) have magnifica- 
tions of about 1.7. Bottom right: Seven-image geometry of source 4, which 
is in a butterfly cusp. The highest magnification image is at (0, —2.9) with 
a magnification of about 7; the next highest are at (±2.9,-1.6) with a mag- 
nification of about 5.6. 
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Figure 9. Histograms of the magnification of the brightest image (left col- 
umn) and the maximum observable time delay between images (right col- 
umn) for the five-image, core triplet and disc triplet configurations for the 
Milky Way model. The effect of removing the halo is indicated by the dotted 
histograms. 



Table 2. Cross sections for the three, five and seven image geometries for 
the Milky Model, the sub-maximum disc, and the Milky Way model with 
halo removed 





Milky Way 


Sub-Maximum Disc 


Halo Removed 


7 image 


0.03 


0.008 


0.02 


5 image 


0.52 


0.23 


0.53 


Total 3 image 


1.57 


1.03 


1.61 


Core triplet 


0.89 


0.85 


0.82 


Disc triplet 


0.68 


0.19 


0.79 



with <E>b, 3>d, and $h the potentials of the bulge, disc and halo 
respectively: 



-GM h 
f + r 



-GM d 



R 2 + (A + Vz 2 + B 2 Y 



(27) 



GM C 



1 , f 2 . r c _i / f \ 



The halo potential corresponds to a cored isothermal sphere with 
density 

Pc 



p(r) 



1 + f 2 /r 2 



(28) 



where M c = 4,-Kp c r^. A typical set of vales for the parameters are 
M b = 3.4 x 1O 1O M , r = 0.7kpc, M d = 1O U M , A = 
6.5 kpc, B = 0.26 kpc, M c = 5 x 1O 1O M , r c = 6.0 kpc. The 
rotation curve of this model is shown in Figure[5] The local circular 
speed is ~ 220 kms -1 and the model has a flattish rotation curve 
out to ~ 50 kpc. 

Let us consider lensing by an edge-on galaxy with these com- 
ponents. (The effect of the disc component decreases rapidly as the 
inclination angle departs from zero). Choosing the scale £o = ro — 
0.7 kpc, we readily obtain 

a » 9.3 , b w 0.37 , — » 2.94. (29) 

m b 

As for the cored isothermal sphere, its dimensionless projected den- 
sity k{t) = a(ror)/a CI is 

-1/2 



: Ph \l + AVI 



where the dimensionless and A are 

TTTc \ — — 

& cr ' c 



ph = Pc 



(30) 



(3D 
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The bending angle is 

r 

Qh = 2ph: 



1 + 



(32) 



We include a halo for completeness, although we will show later 
that it has a small effect on the lensing properties of the model. 



4.2 Typical Image Configurations 

Consider, for concreteness, a galaxy like the Milky Way at a red- 
shift of 0.4 lensing a quasar at a redshift of 1.5. In an Einstein- 
de Sitter universe with Hubble constant of 50 kms _1 Mpc _1 , the 
redshift-distance relation gives cr cr w 2.17 x 10 9 M kpc -2 , and 
the dimensionless mass/density parameters are rab ~ 10.2, md ~ 
30 , ph ~ 0.16 . (Figure[6]shows the dependence of a CI on source 
and lens redshifts in an Einstein-de Sitter universe. As distances 
change, the dimensionless mass/density parameters scale accord- 
ingly). The critical curves and caustics of this lens are shown in 
Figure |7] There are 1, 3 or 5 images for sources in the different 
main regions. Notice that there is a tiny seven-image region within 
the butterfly cusp, shown enlarged in the rightmost panel. The for- 
mation of butterfly and swallowtail cusps in some lens models has 
been noted b efore - for example, in the study of disc-like and boxy 
ellipticals bv lEvans & Witt feooih. an d in the study of the effects 
of halo substructure by feradac et al] J2004h . It is interesting that 
a straightforward disc and bulge system by itself can give rise to 
such h igher order imaging, a result also found by IWang & Turner! 
dl997l) . 

For a fixed source position (£, r? ), the images are at sta tionary 
points of the Fermat surface (see e.g. lSchneider et al] JT992)) 



(33) 



where ip(x,y) is the deflection potential of the three-component 
(bulge, disc, halo) model: 



ip{x,y) = tpb(x,y) +ip d (x,y) +iph(x,y) , 



(34) 



and 

ifa> 

4>h 



-m b (log— +2x(r)), 



2 m d log 



+ (a + y/lP+y*) 1 



(35) 



= 2 /9h A~ 2 [ V / l + A 2 r 2 - log(l + v 7 ! + A 2 r 2 ) 



Here, ipt, and ifa are obtained from the relation 



2 / r 

10 



'/Sj(r')log(— ) dr' 



(36) 



In Figure [8] are plotted Fermat surfaces for four characteris- 
tic source positions. The locations of the images are shown. The 
multiple-image configurations are different to those of spherical or 
elliptical lenses. In all four cases, images near the plane of the disc, 
typical of such lenses, are demagnified. The threefold image con- 
figurations may be split according to the whether the images are 
offset to one side of the galaxy centre ('disc triplets') or whether 
one image occurs close to the galaxy centre and the other two on 
either side ('core triplets'). In fivefold or sevenfold imaging, there 
is one central image whilst the remaining images lie on a broken 
Einstein ring. 



4.3 Magnifications and Time Delays 

Although interesting, the caustics are not directly observable. In 
addition to the positions of images, the observables of any strong 
lens may include the time delays and the ratios of image fluxes. The 
time delay is related to the Fermat potential via 



c-2 
SO 



c D\D\ 



-(1 + zy)4>^ n) {x,y) 



(37) 



where z\ is the lens redshift, here taken by 0.4. To investigate the 
distributions of these quantities, we generate 10 3 source positions 
randomly within the caustics for each of the five-image and three- 
image cases and record the image properties. The magnifications 
are ordered fix, fi%, ■ ■ ■ in increasing order. The means of the least 
highly magnified image, the second least highly magnified image, 
and so on, are recorded in Table Q] whilst histograms of the dis- 
tributions are shown in Figure [9] Similarly, the time delays are or- 
dered 1 1 , t% , ■ ■ ■ in increasing order. For the 5 image and core triplet 
cases, the most delayed image is usually the unobservable, highly 
demagnified central image, which corresponds to ts or tz respec- 
tively. The histograms show the maximum observable time delay, 
which is the largest time delay between the remaining images after 
the central image has been excluded. 

Typically, the largest time delays occur for the core triplet con- 
figurations and next for the five-image systems. For lenses compa- 
rable to the Milky Way, it is of the order of month on average. 
However, the three images in disc triplets have very similar arrival 
times and so the time delays are much smaller, about a day on av- 
erage. The images of highest magnification typically occur in five- 
image systems, for which the average total magnification is 4.7. 
Disc triplets are more highly magnified than core triplets, the av- 
erage total magnifications being 4.2 and 2.1 respectively. Notice 
from Table Q] that disc triplets typically consist of two images of 
very similar brightness and one faint image. By contrast, the av- 
erage magnifications of the three images of core triplets are more 
disparate. 



5 ASTROPHYSICAL APPLICATIONS 

In this section, we consider changes to our basic Milky Way model 
to investigate the lensing characteristics of sub-maximum discs 
(§5.1), the importance of the dark halo (§5.2) and the cross sections 
to multiple lensing (§5.3). 



5.1 Spiral Galaxies and Maximum Discs 

Here, we construct a spiral galaxy model with a sub-maximum disc. 
It has the same bulge as our Milky Way model, but the mass of 
the disc is halved. To maintain the amplitude of the rotation curve 
as ~ 220 kms -1 between 15 and 30 kpc, the core radius of the 
dark halo is adjusted to 4.5 kpc. In the sub-maximum model, the 
contributions of the disc, bulge and halo to the total circular speed 
217 kms -1 at the Sun are 110, 123 and 141 kms -1 , respectively. 
In other words, the dark halo is now the single largest contributor 
to the rotational support at the Solar radius. The rotation curve of 
the model is shown in Fig.|5]as a dotted curve. 

Given the caustics, numerical integration gives the areas of the 
various regions of multiple-imaging. It is traditional to normalize 
the cross sections with that of a singular isothermal sphere with 
circular velocit y 220 km s" 1 . This cr oss section is 7r(£ s j s /£o) 2 f7o j 
where (see e.g., ISchneideret all 19921) 
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can roughly estimate its effect on the cross-sections using the mean 
magnifications listed in Table [JJ Suppose, for example, the lumi- 
nosity function of the sources has a slope of -2. Then, the number 
of five image systems for the sub-maximum disc actually becomes 
similar to that for the Milky Way case. Similarly, the cross-sections 
of the core triplets and disc triplets in the sub-maximum case must 
be boosted by factors of ~ 1.5 and ~ 2.1 relative to the Milky Way 
case to incorporate the effects of magnification bias. Even so, core 
triplets will still predominate strongly over disc triplets. 




12 20 28 36 44 52 

'2 - 'i (days) 



0.28 0.42 0.56 0.70 c 

h ~ 'i (days) 



Figure 10. As Fig. [9] but for a sub-maximum disc spiral galaxy. 



(38) 



and the velocity dispersion i s related to the circular ve locity v c 
b y v c = V2cr v (see e.g., iBinnev & Tremaind Il987h . Just as 
in Keeton & Kochanekl ( fl998h . the magnification bias is neglected 
in the computation of the cross sections. 

The cross sections for the multiple image geometries of the 
sub-maximum disc, normalized to that of an isothermal sphere, are 
given in Table[2] They can be directly compared to the same quanti- 
ties for the Milky Way, which are also listed. For the sub-maximum 
disc, the cross section of disc triplet imaging decreases by a factor 
of ~ 1/3, five-imaging decreases by a factor of ~ 1/2, whilst the 
cross section of core triplets remains largely unchanged. This is 
because the astroid caustic shrinks dramatically in size, whilst the 
radial caustic undergoes a more modest change. In the absence of 
selection effects, a clear-cut difference between maximum and sub- 
maximum discs is that the former gives rise to core and disc triplets 
in roughly equal numbers, while the latter gives rise to mainly core 
triplets. 

Table [JJ shows the magnifications of all the image configura- 
tions, whilst Fig [TO] gives histograms of maximum magnifications 
and time delays. Compared to the Milky Way, the total magnifi- 
cations of all the image geometries has increased. This is because 
the size of the caustics has reduced, and so the source is typically 
close to a caustic when multiple imaging occurs. Although multiple 
imaging is less likely (the cross sections are smaller), the images 
are more highly magnified. Also, the time delays are smaller, as the 
total projected lensing mass is reduced. 

Thus far, we have ignored the magnification bias. However, we 



5.2 The Role of the Dark Halo 

Here, we take the Milky Way model and remove the dark halo. The 
resulting rotation curve is shown in Fig[5]as a dashed line. For the 
case of the Milky Way, the dark halo has an almost negligible effect 
on the size and shape of the caustics. The removal of the halo causes 
changes of ~ 10% in the cross sections of the disc and core triplets, 
as recorded in Table [2] However, the cross sections for total three- 
imaging (i.e., both core and disc) and five imaging are virtually 
unchanged. In fact, the cross section for multiple imaging is very 
slightly increased with the removal of the halo! This is because the 
astroid caustic increases in size slightly more than the radial caustic 
shrinks. 

The effects of the dark halo on the image magnifications and 
the time delays are shown as the dotted histograms in Figure [9] 
whilst the averages are reported in Table [2] Although removal of 
the halo causes a slight diminuition in the time delays, the shape 
of the time delay histograms are virtually unaffected. There are, 
however, some small but noticeable changes in the histograms of 
the magnifications. The total magnification diminishes by ~ 30% 
on removing the halo. 

In our Milky Way model, the contributions of the disc, bulge 
and halo to the total circular speed 225 kms -1 at the Sun are 
155, 123 and 106 kms -1 , respectively. In other words, the dark 
halo makes a modest contribution to the rotation curve and hence 
the mass budget within the Solar circle. The typical scales probed 
by strong lensing are only of the order of a few kpc. Therefore, 
strong lensing by galaxies such as the Milky Way does not probe 
the dark matter distribution in the halo very effectively. We con- 
clude that a Hernquist bulge and Miyamoto-Nagai disc by them- 
selves provide a realistic and analytically tractable lensing model 
for galaxies with maximum discs like the Milky Way. 

5.3 Cross Sections for Multiple-imaging by Spiral Galaxies 

Figures[T71andll2lshow cross sections of multiple-imaging varying 
with the masses of disc and bulge, and with the shape parameters a 
and b of the two-dimensional disc model. The cross sections depend 
on the distances between lens, source and observer through a cr , 
but changing a cr simply corresponds to rescaling the dimensionless 
mass/density parameters m^, m^, p^, so cross sections for a CI 7^ 
2.17 x 1O 9 M0 kpc -2 can also be read off from Figure[6] noting 
also that the singular isothermal sphere cross section, from 138t . 
goes as a~ 2 . In all the panels, the position of the Milky Way galaxy 
is shown by a solid circle. 

From Figure[TT] we can read off the cross sections for a galaxy 
like the Milky Way (see also Table |2j- The cross sections show 
that the core triplets, disc triplets and fivefold imaging are roughly 
equally likely. Seven imaging configurations are roughly a thou- 
sandfold times less likely. Increasing the mass of the bulge Mb at 
fixed disc mass Afd causes the cross section for disc triplets to di- 
minish and that for core triplets to increase. Increasing the mass of 
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Figure 11. Multiple-imaging cross sections, as a fraction of the singular 
isothermal sphere cross section, varying with bulge and disc masses. Masses 
are in solar masses (the Milky Way values are in the centre of each plot). 
All other parameters are set to Milky Way values. Disc to bulge mass ratios 
are constant along diagonals. Top left: Disc triplet imaging cross section. 
Top right: Core triplet imaging cross section. Bottom left: Five-imaging 
cross section. Bottom right: Seven-imaging cross section. The butterfly cusp 
cannot be made arbitrarily large by making masses larger. 



the disc Md at fixed bulge mass Mb causes the cross sections for 
fivefold and sevenfold to increase. 

To interpret Figure[T2] we recall that that the parameter a con- 
trols the ellipticity (a = is circularly symmetric), while b controls 
the sharpness of the disc profile, 6^0 corresponds to a razor thin 
disc). Increasing b at fixed a causes the cross sections of disc triplets 
and quintuplets to fall sharply, whilst leaving the cross section for 
core triplets largely unchanged. 



6 CONCLUSIONS 

We have analyzed the strong lensing properties of a realistic model 
of the Milky Way galaxy with a disc, bulge and halo, combined 
to produce a nearly flat rotation curve. All three components - the 
spherical Hernquist bulge, the Miyamoto-Nagai disc and the cored 
isothermal halo - have analytic deflection angles. There is strong 
evidence that the Milky Way galaxy possesses a maximum disc, in 
the sense that the contributions to the rotation curve by the disc and 
bulge at the Solar radius dominate over the contribution of the halo. 
The consequence of this is that the halo has a small effect on the 
strong lensing properties, including the magnifications and the time 
delays. 

As a strong lens, the Milky Way galaxy exhibits four different 
kinds of multiple imaging geometries. They are (i) three images 
on one side of the galaxy centre ('disc triplets'), (ii) three images 
with one close to the center ('core triplets'), (iii) five images and 
(iv) seven images. For spiral galaxies, the lensing cross sections are 
dominated by edge-on models, as this is where the disc makes its 



Figure 12. Multiple-imaging cross sections varying with the disc parame- 
ters a and b. All other parameters are set to Milky Way values. The four 
panels correspond to the same cross sections as in Figure [T71 



presence known most effectively dKeeton & Koch anek 1998j). The 
edge-on case is also the most important from the point of view of 
applications. In this instance, the cross sections show that the core 
imaging, disc imaging and fivefold imaging are roughly equally 
likely. The disc triplet cross section is also sensitive to the thickness 
of the disc and the mass of the bulge. 

Disc triplets are a characteristic of gravitational lensing by 
spiral galaxies with maximum discs. They occur when the astroid 
caustic pierces the radial caustic. They consist of three images 
straddling one side of the lensing galaxy. The two outermost images 
are usually of comparable magnification, whilst the intervening im- 
age is much fainter. All three images have very similar arrival times 
and so the time delays are small, less than a day on average for our 
representative model. For comparison, the largest time delay be- 
tween the observable images in core triplet or fivefold imaging is 
of the order of month on average. 

Spiral galaxy lenses are rare compared to early-type galaxy 
lenses. The importance of gravitational lensing by nearly edge-on 
spirals is that it directly probes the shape of the matter distribu- 
tion and therefore can break the degeneracy between flat disc and 
round halo. Face-on spirals are ineffective in discrimination, as the 
projected matter distribution from the disc is round. Even though 
only a few examples of nearly edge-on spiral galaxy lenses are 
known, this has already pr ovided some important results. For ex- 
ample, iMaller et al.l ( feOOfj) has concluded that the spiral lens of 
B 1600+434 does not have a maximum disc from detailed mod- 
elling of the positions and flux ratios of the two visible images. 

The lensing properties of spiral galaxies with sub-maximum 
discs differ from maximum discs. The astroid caustic shrinks sig- 
nificantly in size, whilst the radial caustic shrinks more modestly. 
Therefore, the cross section of disc triplets is substantially reduced, 
whilst the cross section of core triplets remains roughly unchanged. 
The overall effect is to reduce the cross section to threefold (disc 
and core triplets) imaging by ~ 2/3 and to fivefold imaging by 
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1/2. The total magnification of the configurations is on average in- 
creased, but the time delays are decreased. 

If maximum discs are typical, the lensing cross sections sug- 
gest that core triplets predominate only slightly over disc triplets. 
However, if sub-maximum discs are typical, then core triplets are 
~ 9/2 times more likely than disc triplets. In fact, except for 
APM08279+5255, no disc triplet configurations are known. One 
possible interpretation of this scarcity is that the Milky Way is 
atypical and that spiral galaxies are usually sub-maximum. This 
conclusion derives some support from the detailed modelling of 
iMaller et all j2000h in the case of B 1600+434. 

It is also interesting that sevenfold imaging can occu r for disc 
galaxies (c.f. IWang & Turned [l997l : lEvans & Witdl200lh . even if 
the cross section is lower by a factor ~ 10 -3 compared to the cross 
sections for threefold or fivefold imaging. Our cross sections do not 
incorporate the effects of magnific ation bias, which will boost th e 
likelihood by a significant amount dTurner. Ostriker. & Gottll984l) . 
Note that the existence of higher order cusps increases the magni- 
fication of the quintuple image configurations as well. This is be- 
cause the total magnification of the 5 images is increased by the 
appearance of the butterfly cusps inside the 5 image caustic. This 
fact has been neglected in all calculations of the numbers of ex- 
pected 3 or 5 lens systems in the literature. 
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